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Abstract 

Diffractive production of J/\I/ particles by virtual photons on a pro- 
ton target is studied with a view towards understanding two important 
corrections to the leading order result. First, the effect of Fermi motion 
of the heavy quarks is studied by performing a systematic expansion in 
the relative velocity, and a simple correction factor is derived. This is 
considerably less than estimated previously. Second, since the kinematics 
necessarily requires that non-zero momentum be transferred to the pro- 
ton, off-forward gluon distributions are probed by the scattering process. 
To estimate the importance of the off-forwardness, we compute, in leading 
order perturbation theory, the extent of deviation from the usual forward 
gluon distribution in a quark. 



1 INTRODUCTION 



There is considerable interest in diffractive electroproduction of J/^f mesons off 
protons at high energies because this process is important for studying the gluon 
density in a proton at low values of Bjorken x. This interest stems from the sim- 
plicity of the leading order QCD expression for longitudinally polarized photons 
which was first derived by Ryskin |2[. 

^( 7 * + p - J/* + p) = 16 ^ Tee *s(Q 2 )MZ, Q 2 )} 2 (l) 

where, Q 2 = \(Q 2 + M 2 ), M is the J/^> mass, and T ee is the decay width into 
leptons. The above equation was derived under the assumptions of s 3> Q 2 3> 
M 2 3> t, and that Fermi motion of the quarks in the meson can be entirely 
neglected. It was further supposed that the gluon density appearing in eq.|l| 
is that which would be measured in some inclusive hard process, i.e. that it 
corresponds to the matrix element of gluon operators between states of equal 
momentum. 

In this paper we shall examine the effect of relaxing two assumptions which 
go into eq.|l|. The first is to take into account the correction arising from the Fermi 
motion of the cc pair. In the work of Brodsky et al[|J this motion is contained 
in the vector meson light-cone wavefunction ip v (k±,x), a quantity which is in 
principle calculable from lattice QCD but whose presently unknown form is an 
important source of ambiguity. For example, Frankfurt et al|§ conclude that 
wavefunction effects suppress J/\l/ production by a factor of 3 or more. However 
Ryskin et al[Q] estimate a suppression factor of 0.4 < F 2 < 0.6. The detailed 
shape of the wavefunction appears to be an important source of the difference. 

The method of treating the diffractive process, as well as Fermi motion 
corrections, used in this paper differs from previously used methods in an essential 
way. Rather than work in the infinite momentum frame and in the A + = 
gauge, we shall choose the rest frame of the J/\l/ and the Coulomb gauge for 
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the soft gluons in the meson wavefunction. This is the natural choice for heavy- 
quark systems because one can then use systematic procedures, such as Non- 
Relativistic Quantum Chromodynamics|| (NRQCD) or the method developed 
in refs. [7|, ||, in order to evaluate quarkonium observables of interest to any 
desired level of accuracy. However, for the gluons in the fast moving proton we 
shall continue to use the A + = gauge because this is the natural gauge to use 
for parton distributions. The two kinds of gluons have very different momenta 
and hence are effectively distinguishable, justifying the use of two different gauges 
in two different parts of the same Feynman diagram. It turns out that a gauge- 
invariant correction factor, derived in this paper, multiplies eq.l, 

V 9M 2 J 

The second derivative of the wavefunction is understood to be evaluated at the 
origin. It is a non-perturbative quantity whose value has to be inferred from 
some other quarkonium processes, such as decays or production, involving large 
momentum transfer. In previous work|| [7| its value was estimated, 

V 2 

« -0.07 

M 2 

The correction factor due to Fermi motion is therefore around 0.96, a value con- 
siderably below the other estimates]^, [3J. Hence, the ambiguity in extracting 
the normalization of the gluon distribution may be under better control than 
anticipated so far. 

The second issue to be considered in this paper is the gluon distribution 



which appears in eq.l. Recently JiflTOj identified certain twist-two "off- forward" 
quark distributions inside the proton which, when measured, will reveal the or- 
bital angular momentum content of the proton. Subsequently Radyushkin|TT|| 
extended the discussion to the off-forward or "asymmetric" gluon distribution 
in the proton and pointed out that diffractive vector meson electroproduction 
necessarily measures this quantity. Here we have examined this issue further by 



considering gluon radiation from a quark and explicitly computed the off- forward 
gluon distribution in a quark to leading order in a s . This enables an estimate to 
be made of the extent to which the gluon distribution measured in J /if! diffractive 
production differs from that which would be measured in some inclusive process 
like 7 + P -> J /iff + X. 



2 FERMI MOTION 

2.1 Kinematics 

We consider a massless proton target and the t — limit. Define two null vectors 
and with p 2 = n 2 = and p ■ n = 1, 

p" = A(l, 0,0,1), „" = --L(l, 0,0,-1) (2) 

pt 1 is also the proton's momentum. Although we shall not need to do so explicitly, 
A can be adjusted to bring the produced J /iff to rest. The kinematic region of 
interest is considered to be s ^> Q 2 ^> M 2 . With the definitions, 

(3) 





a 2 

y 2 


= -Q 2 , 


are, 






<f 








= p^ + A^ 




A 11 








= 





(4) 

The polarization vectors of the longitudinally polarized photon and J /iff are, 
respectively, 



These obey e L .e L = 1, e L .q L = 0, and E L .E L = -1, = with K 2 = M 2 . 

We have kept only leading terms and set Aj_ ~ 0. 

2.2 Diagrams 

The leading order contribution to J/\I/ diffractive production is given by the sum 
of the diagrams shown in Fig. 2, to which must be added the contribution of two 
other diagrams that give the same numerical values because of time-reversal sym- 
metry. Consider, by way of example, the first of these which has the expression, 

Ai = / (^i ^- 4 Tr[S, v (k,A)Hr(k,£)M(£)), 
HF(k,£) = e Q g 2 YS F (k + q-K/2 + £)YS F (q-K/2 + £)/(q). (6) 

The perturbative part H^ u (k,€) is different for the other diagram but the other 
factors in eq|] remain unchanged. We have not indicated colour explicitly in the 
above; its inclusion will amount to a simple factor which will be inserted at the 
end of the calculation. The non-perturbative information of the vector meson is 
contained in the Bethe-Salpeter wavefunction M{£), 

M{£)= J d 4 xe i£ - x (K,E\T[tP(x/2)^(-x/2)]\0). (7) 

In the above, £^ and x M are, respectively, the relative momentum and relative 
distance of the cc pair. The non-perturbative information of the gluons in the 
proton is contained in S^, 

S^(k,A) = J d 4 xe i{k+1 ^ A) - x (P'\T[A^-x/2)A u (x/2)]\P) (8) 

While the diagrams in Fig. 2 contain the leading order contribution to the cross- 
section, they also contain parts which are next to leading order (NLO). The sense 
in which these are to be understood as "higher order" will be made precise later. 
Other diagrams will have to be included (see Fig. 3) for a complete calculation at 
the NLO level. 
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2.3 Expansion 



The diffractive process considered here has two large scales, Q 2 M 2 3> Aq CD . 
Since a cc system is close to being a non-relativistic coulombic bound-state, it 
allows for an expansion in powers of the heavy quark relative velocity. Hence it 

is useful to expand the inner integral in eq.^, 

r rJ 4 £ 00 

n W = / 7^^(M)M(£) = £ njr 

00 1 <9 9 

= Err^r---^r^l^M-- (9) 



where, 



„ n\d&* ()('»* 



M ai~*« = 1 -0^£ ai ■ ■ ■ £ an M(£) 



= id ai ■ ■ ■ id an (K, E\T[if>(x/2)ip(-x/2)]\0) \ x=0 (10) 

The set of constants M ai "' an provide a description equivalent to that of the 
original BS wavefunction in eq.|7|. The expansion eq.|9] is useful because the quarks 
are nearly on mass-shell: (~if + £) 2 ~ m 2 implies that all components of are 
small relative to the quark mass m in the meson's rest frame and, in particular, 
£ ■ n ~ (m/Q) 2 . In the large Q 2 limit this implies considerable simplification, 
giving a limit approximately independent of £, 

! „ K \ 

(k + q - K/2 + £) 2 - m 2 + is Q 2 k ■ n - £ + is v ; 

1 2 
{q - K/2 + £) 2 -m 2 + ie ~ 'Q 2 ^ 

Inclusion of Fermi motion requires that we keep a sufficient number of derivatives 
w.r.t £ in eq.|9|. These may be computed using the simple Ward identity, 

JL Sf = -S Fl a S F) (13) 

and the Q 2 — > oo limit should be taken after performing the trace algebra. Stated 
in words, a differentiation of either propagator in eq.^ w.r.t £ splits that propa- 
gator into two. Since we shall work upto 0(v 2 ), only two derivatives of H^{k 1 £) 
are needed. 



5 



2.4 Gauge Invariance 



It is obvious from the occurence of the ordinary derivatives in eq.[TO, or the form 
of the BS wavefunction eq.0, that gauge invariance has been violated. In earlier 
work on quarkonium processes 0, [7|, || we have encountered an identical situation 
- the diagrams in Fig. 2 yield expressions which are not gauge invariant to 0(v 2 ) 
and one needs to consider additional diagrams, which are higher order in a s . 
These are shown in Fig. 3. The gluon fields indicated in these diagrams combine 
with the ordinary derivatives to yield covariant derivatives, d a —>■ D a , thereby 
restoring gauge invariance. In the Coulomb gauge, the contribution of explicit 
gluons is 0(v 3 ) and so the reduction of the Bethe-Salpeter equation performed 
in ref.0 without explicit gluons is adequate upto 0(v 2 ). We therefore arrive at 
the following gauge-invariant matrix elements: 

6 M 2 \ M) 

In the above, <fi and V 2 are the non- relativistic wavefunction and its second 
derivative evaluated at zero separation. Inclusion of V 2 amounts to taking the 
first step towards inclusion of Fermi motion. 

2.5 Traces 

All the ingredients are now in place for calculating the trace of the quark loops. 
Because we shall need only the leading twist piece, symmetric in /i and u, it will 
be sufficient to calculate, 

n n = TrJ2 K = (-9^ + V»n v + p v n^)Tr[^} (15) 

i=l,2 



for n — 0, 1,2 (n is the order of differentiation w.r.t £) and then keep only the 
leading order term in 0(1/Q). We record below the results of the calculation 
listing, for clarity, the relative contribution of only those diagrams which give a 
non-zero contribution, 

* - -^* + * (i + £R) 

ft! = 0(1/Q 3 ) 

" 2 ~ M^Q 2( 3 + 3 + 3~3 ) MV +0(Vg) (16) 

The factor of 2 multiplying the brackets in the above equations comes from the 
diagrams which are permutations of the ones shown. The sum over all diagrams 
is, 

--^)( 1 + ^) +0(1/e3) ' (17) 

Note that this leading order contribution is in fact independent of the gluon 
momentum k in the Q 2 — > oo limit. The term in the brackets represents the 
correction due to the Fermi-motion of the heavy quarks and its square is precisely 
the factor which modifies eq.l. 



3 GLUON DISTRIBUTION 

3.1 Asymmetric Distribution 

Let us now return to the amplitude for diffractive scattering, a typical contribu- 
tion to which is given by eq.|| The photon and proton both move along the z 
direction, and the gluons in the proton have limited k\ and k 2 . This means that 
one can perform a systematic collinear expansion in these quantities just as in 



the treatment of deep-inelastic scattering[12 



H^(k,£) = H^(k + ,£) + (k — k + ) a d a H^(k,£) \ k=k+ + 
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Keeping only the first, leading twist, term gives in the A + = gauge, 

r r d\ \ \ 

= jdyj ^e^-^\P'\A^n)A u qn)]\P) H^(y,£) (19) 

In the above we have set x~ = Xn~ and k + = yp + . The time ordering operation 
becomes irrelevant on the light- cone. 

The inner integral will now be analyzed following the discussion given by 



Radyushkin[ll]. Define the "asymmetric distribution function", F^(X), as below, 



(P'\n-G + \-±n)n-G+\±n)]\P) 



= u(p')$u{p) dX {e^ x -^ + e- iX ( x -W} F^X) (20) 
o 

A sum over transverse components (i = 1,2) is implied. The proton spinor 



product is u(p') fiu(p) = 2 a/1 — £, with the initial and final protons having the 
same helicity and p 1 = (1 — £)p. Making a Fourier transformation yields, 

F&) + m -V) = 7~^= / d ^y-^\P^n-G^\-^n)n-G^\^n)]\P) 

(21) 

It is instructive to insert a complete set of states, 

Ft(y)+m-v) = -77=4 £ -l+*) + 8(Z-V-l + x)} (P'\A*\k)(k\A*\P) 
4 vi-4 k 

(22) 

Here x = k ■ n with < x < 1 is the momentum fraction carried by the interme- 
diate state. Comparing with the usual (diagonal) gluon distribution function for 
£ = it immediately follows that, 

F e=o(y) = \ya(y), 

9(y) = y'£6(y-l + z){P\A i \k)(k\A i \P) (23) 
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We shall now relate the matrix element in eq.|T^ to F^(y). Inverting the 
relation G +l = d + A 1 gives, 



A^Xn) =rT da G +i (Xn + an) 



Inserting the above into eq.[19] and using the definition of F^(y) in eq.|20 

dX 



J^ e iHv-^ {P )A\-^n)A^n)\P) 



+ 



y(€-y-ie) - y - ie)(y - ie) 
The imaginary part of the above for y > is 



and hence, 



(24) 



(25) 



Im I ^-0^Tr[S^(k,A)H^(k,£)M(£)} = . (26) 



3.2 Cross-section 



All the ingredients are now in place for calculating the cross-section for the diffrac- 
tive process under consideration, 

1 



da 
~dt 



16ns 2 



A 



1 



(27) 



The factor comes from summing over colours, and Q is the quantity calcu- 



lated in the previous section, eq.[L7| , from the expansion of the heavy quark loop 
integral. Defining T to be the leading order decay width into lepton pairs, 



2 „,2 



M 2 ' 



yields the following important result 
da 16vr 3 Mr 



dt 



3Q 6 



-a s (Q 



1 + 



8V 2 \ 



(2* 



(29) 



Making the approximate identification, 



4^1-^(0^(0, (30) 
and setting the last factor to unity reproduces eq.l once again. This identification 



was motivated by eq.[23| but the exact relation between and g(0 is far from 

clear. 

3.3 Perturbative gluon distribution 

and g(£) can be known only if the non- perturbative structure of the proton 
state is known. However, it would be highly desirable to have at least some partial 
knowledge of their structure. To this end, consider the following simple solvable 
problem: imagine that the target proton is replaced by a single quark which can 
radiate a gluon. Its light-cone wavefunction can be computed order by order in 
perturbation theory, and the leading order matrix element is, 

(?s\A^ks) = g^±^u( P , s'Xl\)u(ps)e* i (l\), (31) 

where I, A are the momenta and transverse polarizations of the emitted gluon. 
Using, 

£^AK*(/A) = -^+ , + , (32) 
and summing over i = 1,2 gives, 

(pr\A*\k)(k\A*\P) = g *-^=±. 1 + * 2 ZA (33) 
Inserting this into eq.E2 yields, 



2 y{£-y) [ d2k 



Wy -l + g)+ ^- y -l + x)} (1 _ 1 + ( f_-/_ 

<*. f 1 , (l-y) 2 + (l-e + Z/) 2 ) [dkj 

4 7 r\ 1+ 2(1-0 \J kl [6V 
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The last integral is both infrared and ultraviolet divergent. It is regulated by 
inserting a low momentum scale cutoff [i ~ O(Aqcd) and a high momentum 
cutoff k± ~ 0(Q). Multiplying by the colour factor Cp = |, we arrive at the 
perturbative asymmetric gluon distribution inside a quark, 

*»-£{ 1+ %^?M (35) 

Note that, 

Vi-e^(0 = ^ i v /rr ^{ 1 + ( 1 -o}iog^ 

but that the usual perturbative symmetric distribution, which can also be ob- 



tained by first putting £ = in eq.[35| and then setting y — £ is, 

fcKfl = ^(i-e+^ a V (37) 



3vr V "2 

Comparison of the last two formulae gives an estimate of the extent to which the 
asymmetric distribution departs from the symmetric one as £ becomes larger. 

Finally, we remark that there exists some confusion in the literature about 
various factors of four. First, it is claimed in the work of Brodsky et al|4[] that 
the cross-section displayed in eq.l must be multiplied by \. We do not find this 
to be the case; the result of Ryskinfl], || appears to be correct. A second point 
of confusion, unrelated to the first, is contained in the work of Radyushkin||l 1|| 
who states that the usual gluon distribution, is approximately related to the 



asymmetric distribution F^(X) by the relation y/1 — £Ff (£) — > £<?(£), i.e. without 
the factor of 4 contained in eq.[23] of this paper. Again, we do not agree with this. 
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Figure Captions 



1. Definition of kinematic variables for J/ty diffractive production off a proton 
target by a virtual photon. 



2. Diagrams which give non-zero contribution at order Q~ l and v°. The relative 
weight at this order a:b is as 1:2. Two other diagrams, which are numerically 
equal by time-reveral invariance, are not shown. The complete expression is given 
in eq.[U| 



3. Diagrams which give non-zero contribution at order Q~ l and v 2 . The crosses 
denote connection to external gluons originating from the proton. The relative 
weight at this order a:b:c:d:e:f is as -l:l:-2:2:2:-4. Note that each internal gluon 
zero- momentum gluon line, in the Coulomb gauge, is actually just a differenti- 
ation of the quark propagator. Six other diagrams, which are numerically equal 
by time-reversal invariance, are not shown. The complete expression is given in 
eq.|l6 
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